Staggered currents in the mixed state 
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The current pattern in the mixed state of high-Tc superconductors is studied in the U(l) mean 
field theory of the t-J model. Our findings are the following. 1) In the absence of antiferromagnetism 
a robust staggered current pattern exists in the core of vortices if the doping is not too high. 2) At 
a fixed doping and with increasing magnetic field, the size of the staggered current core expands, 
and eventually percolates. 3) The polarity of the staggered current is pinned by the direction of 
the magnetic field. 4) Vortex cores locally modify the hole density - in a staggered (non-staggered) 
core, the excess charge is slightly negative (positive). 5) Gutzwiller projection does not wcish out the 
staggered current. Finally we present two experimental predictions concerning neutron scattering 
and STM spectra that capture the signature of the staggered current induced by the vortices. 

PACS numbers: 74.25.Jb, 79.60.-i, 71.27.+a 



Recently, there has been a renewal of interest in the so- 
called staggered flux phase of the t-J model 0-0] where 
a circulating current produces a staggered orbital mag- 
netic moment. According to the mean-field theory, 
the staggered flux phase, although never stable, is a close 
competitor of d-wave superconductivity when antiferro- 
magnetism is absent. This situation can be aptly de- 
scribed by a Ginzburg-Landau (GL) theory with two 
competing order parameters. Q A serious concern is 
whether the staggered flux phase actually exists in the 
high-Tc phase diagram. |^||] 

Even if the staggered flux phase does not exist, GL 
theory would predict that in the core of a vortex, where 
the superconducting order is suppressed, the staggered 
flux order has a chance to appear. This is first pointed 
out forcefully in a recent paper by Lee and Wen jl^ from 
the viewpoint of their SU(2) mean field theory of the t- 
J model. Recently we have shown that the staggered 
current core also exists in the U(l) mean-field solution 
of the same model. Moreover we demonstrated that 
such current pattern survives the Gutzwiller projection 
which removes any double occupation in the mean field 
results. 

The purpose of this paper is to address several impor- 
tant issues concerning the staggered current in the mixed 
state of high-Tc superconductors. One of our main results 
is the prediction for two experiments which are relevant 
to the existence of staggered current in vortex cores. Our 
conclusion is based on the U(l) mean-field solution and 
its Gutzwiller projection. 

Our starting point is the t-J model with Coulomb in- 
teraction, 

^ = E (CL^^..e-^- + h.c.) + JE(S. • S, - ^) 
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In the above Aij = (27re//ic) A • dl is the link phase 
produced by the physical vector potential A. All other 
notations are standard. In the rest of the paper we use 
the parameter set t = V = 3 J. 

In the U(l) slave-boson approach the following replace- 
ments are made in Eq. (|l|): 1) Cj^Cja ^ fj„fjab]bi, 2) 
S, (l/2)/t^fTc/3/j/3, , 3) ^ 1 - bib,, and 4) the 
occupation constraint — * + J2a fiafic^ — 1- 
above hi is the boson and fi^ is the fermion annihilation 
operator respectively. Our calculation is performed at 
zero temperature where the bosons are assumed to have 
condensed. The spin-exchange term is replaced by the 
following mean-field decoupling 
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where Kij and Ay , the mean-field hopping and pair- 
ing amplitudes, are determined self-consistently, together 
with the condensate boson order parameter (bi) and the 
Lagrange multiplier that enforces the occupancy con- 
straint. Due to the internal U(l) symmetry we can al- 
ways choose the gauge so that (bi) is real and positive. It 
is important to stress that the above decoupling excludes 
the spin magnetic moment (S^) as an order parameter, 
and hence will not be able to describe low doping regime 
in which commensurate/incommensurate spin magnetic 
order exists. Aside from the above restriction, no other 
constraint is placed on the mean-field parameters. The 
actual calculation is performed numerically on a x Ny 
lattice under twisted boundary condition. Much of our 
treatment of the vortex lattice is the same as that in Rcf. 

flux 



||l2| . The average magnetic flux density is f — 
quanta, in units of hc/e, per plaquette. 
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The flux density / causes the pairing amphtude Ay to 
wind by 47r around each x Ny unit cell. The question 
concerning whether it is energetically more favorable to 
nucleate two hc/2e vortices or a single hc/e vortex in the 
unit cell has been raised in the literature. [|l0|,^ Our 
result shows unambiguously that it is the hc/2e vortex 
that is favored. More specifically we have checked that 
even when the initial condition corresponds to a single 
hc/e vortex, i.e. two hc/2e vortices on top of each other, 
the final self-consistent solution always exhibits two sepa- 
rated hc/2e vortices. We emphasize, however, that such 
conclusion will be subject to change if the strength of 
the Coulomb potential is modified. The reason for that 
is because the vortex core is charged! 

In Figs. 1(a) and (b) we show two dimensional map of 
the hole density in the unit cell (with two vortices) at 
doping levels x = 10%(a), and 15%(b). The blue region 
marks lower and the red region marks higher hole density 
respectively. Thus the vortex core is negatively charged 
in (a) and positively charged for (b). The charge differ- 
ence between the x = 10% and the x — 15% vortices 
turns out to be related to the presence/ absence of stag- 
gered currents within the vortex core as we will discuss 
below. 
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FIG. 1. Two vortices in a 24 x 24 lattice, (a) Hole dis- 
tribution with X = 10%; (b) Hole distribution at a: = 15%; 
(c) Bond current pattern corresponding to (a); (d) Bond cur- 
rent pattern corresponding to (b). The crosses in (c) and (d) 
highlight the vortex cores. 

In Fig. 1(c) we present the current pattern near the core 
of the lower- left vortex in Fig. 1(a). The variable-length 
arrows represent the direction and magnitude of the bond 
current given by 



-iAi 



(3) 



(Recall that in our gauge the boson condensate amplitude 
{hi) is real and positive.) The pattern clearly indicates 
the existence of a staggered current core. In Fig. 1(d) we 



present the current pattern near the core of the lower- 
left vortex in Fig. 1(b), where the staggered current has 
disappeared. 
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FIG. 2. Bond current patterns at (a) / = 1/576, (b) 1/144, 
(c) 1/64, and (d) 1/36. 

Let us return to Fig. 1(c), and ask what is the correla- 
tion, if any, between the vorticity of the current pattern 
and the direction of the magnetic field. The answer is 
that the circulation around the central plaquette where 
the vortex center reside is always opposite to that of the 
magnetic field. (This is in contrast to the case in 
Fig. 1(d).) Moreover, we intentionally start with an initial 
condition where the circulation of the central plaquette 
is the same as that of the magnetic field, and find in the 
final self-consistent solution that the vortex center moves 
by one lattice spacing, so as to make the circulation of the 
central plaquette and the magnetic field opposite. Thus 
the magnetic field pins the polarity of the staggered cur- 
rent. 

In Figs.2(a)-(d) we illustrate the evolution of the stag- 
gered current pattern found at 10% doping over a wide 
range of magnetic field. The field strength is (a) / = 
1/576, (b) / = 1/144, (c) / = 1/64, and (d) / = 1/36, 
respectively. In physical units the fields considered here 
are very large (/ = 1/1600 roughly corresponds to a mag- 
netic field of 10 Tesla). It is clear that as the vortices get 
closer their staggered current cores overlap and permeate 
the entire lattice, as shown in Figs. 2(c) and (d). ||l^ The 
end result is a state with uniform staggered current. This 
is a magnetic field induced co-existing d-wave supercon- 
ducting and staggered flux state. 

Although the lowest field considered in this paper 
(roughly 30 Tesla) is still very high from the experimental 
standpoint, none of our experimental predictions (as will 
be discussed later) will change qualitatively at a lower 
field. 

Just as in our earlier findings, lO] the staggered cur- 
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rent in the mixed state also survives Gutzwillcr projec- 
tion. Q As an example, Fig. 3 shows the staggered cur- 
rent pattern at x = 6.25% and / = 1/64, before (a) and 
after (b) the Gutzwiller projection. Although the stag- 
gered current is weakened by the projection, it certainly 
does not disappear. 
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FIG. 3. Bond current pattern at x = 6.25% and / = 1/64. 
(a) Mean field result, (b) Result after Gutzwiller projection. 
The crosses mark the locations of the vortex cores. 



quette). In this case the vortex lattice is commensurate 
with one of the bi-partite sublattice of the underlying 
crystal structure. We expect elastic neutron scattering 
will show a similar pattern. 

Seeing the magnetic field induced Bragg peaks near 
(tt, tt) might be also consistent with vortex cores being 
spin-antiferromagnetic. Although it would be dif- 
ficult to distinguish between the staggered orbital mo- 
ments and staggered spin moments from the standpoint 
of neutron scattering, the latter is much less likely from 
theoretical considerations. Due to the presence of the 
external magnetic field and the vortex lattice, all sym- 
metries that ensure the degeneracy between the two po- 
larity of the staggered flux order parameter are broken. 
Thus the staggered flux in the vortex core is an induced, 
rather than spontaneous, order. This is not true for the 
spin antifcrromagnetic order. Indeed, under the experi- 
mental condition, rotation of spins around the magnetic 
field direction remains an unbroken symmetry. As a re- 
sult spin antifcrromagnetic alignment in the vortex core 
requires a spontaneous symmetry breaking, which is not 
possible for a finite (small) system such as the core of a 
vortex. 





FIG. 4. |mkp as a function of k corresponding to Fig. 2(a). 

What are the experimental consequences of the 
staggered-current cores? Figures 1(c) and 2 suggest that 
in a neutron scattering experiment there will be a mag- 
netic field induced Bragg peak near or even at (tTjTt). 
(Whether there is a peak precisely at (tt, tt) depends on 
whether the vortex lattice is entirely contained in one of 
the bi-partite sublattices of the underlying Cu02 plane. 
The peaks around (7r,7r) should be more robust.) To 
make sure that the scattering near (tt, tt) is induced by 
the vortex core, one can 1) study the dependence of the 
intensity and position of the Bragg peaks on the magnetic 
field and 2) correlate the peak position near (tt, tt) with 
that near (0,0), which determines the structure of the 
vortex lattice. ||l5[] In Fig. 4 we present |r7ik|^ , where rrik 
is the Fourier transform of the lattice curl of the current 
(i.e., the directed sum of the bond current around a pla- 
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FIG. 5. (a) Schematic plot of the vortex lattice, and the 
the path along which the density of states is calculated. The 
path begins with one of the vortex cores (crosses), (b) The 
density of states with x = 10% along the path (starting from 
the front), (c) The same as (b) for x = 15%. 

In addition to the scattering experiment staggered vor- 
tex core might also cause detectable difference in the local 
tunneling spectra. In Fig. 5 we present the local tunnel- 
ing density of state as a function of location. The arrows 
in (a) indicate the path (crosses mark the location of 
vortices) along which the local density of states is com- 
puted. Figure 5(b) is the result for x = 10% where there 
is staggered current in the vortex core. Figure 5(c) is 
the result for x — 15%, where staggered current does not 
exist. The most significant difference occur at low ener- 
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gies. While there is a zero-bias peak in the core of the 
15% vortex, there is no such peak in the core of the 10% 
vortex. ||l^ Since the staggered current is known to open 
up a pseudogap the result is not surprising. 

In general spin polarized tunneling experiment can be 
used to differentiate orbital versus spin antiferromag- 
netism. In the presence of the latter, the spin-dependent 
tunneling density of state should show a two-sublattice 
structure. We do not expect such effect for orbital anti- 
ferromagnetism, because a) the Zeeman splitting caused 
by the orbital moment is extremely small, and b) on a 
given site the staggered magnetic field originated from 
four neighboring plaquettes tends to cancel. 

If the signatures of the staggered current core discussed 
above are seen in both the neutron and STM experiments 
it will constitute an extremely strong evidence for the ex- 
istence of staggered current in the vortex core. 

Finally, what do we learn from the existence/non- 
existence of the staggered current in the vortex core? 
From the beginning of high- Tc physics, the t-J model has 
been identified as the model that captures the competing 
local interactions, i.e. charge hopping and spin antifer- 
romagnetic exchange, of the cuprate materials. During 
the last fifteen years important progress has been made 
on numerically simulating this model. Nonetheless ques- 
tions such as whether the ground state is homogeneous, 
and whether it is d-wave superconducting still remain 
controversial. On the other hand the mean-field theory 
of the t-J model has produced very tantalizing results. 
For example, it predicts that d-wave pairing is among 
the most pronounced ordering tendency of the model. 
psf Recently, to a limited extent, mean-field theory has 
also predicted the presence of stripe inhomogeneity. | [Tg| ] 
In our opinion these successes give enough motivation 
to check whether another robust prediction of the mean- 
field theory, the close competition of the staggered flux 
phase with d-wave pairing, is actually correct. The pres- 
ence of staggered current in the core of a vortex is a direct 
consequence of this competition. If such an ordered cur- 
rent pattern is observed, it is reasonable to argue that 
we have understood the rudimentary ordering tendency 
of the model. On the other hand, in the face of negative 
experimental results we should seriously worry about the 
validity of the mean-field theory, or perhaps even about 
our understanding of the relevant local interactions. 
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